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• I review lattice QCD results for hybrid mesons, including a discussion of their 
QO ' hadronic decays. 

O ' 

CO 

. 1. Introduction 

O ■ Quantum Chromodynamics has emerged as the unique theory to describe 

' , hadronic physics. It is formulated in terms of gluonic and quark fields. The 

| only free parameters are the scale of the coupling (usually called Aqcd) 

• and the quark masses defined at some conventional energy scale. 

Where large momentum transfers occur, the effective coupling becomes 

weak and a perturbative treatment is valid: in this domain the theory 
has been tested directly by experiment. However, because the effective 
coupling is weak for these processes that can be described by perturbation 
theory, they are necessarily not the dominant hadronic processes. A typical 
hadronic process will involve small momentum transfers and so has to be 
treated non-perturbatively. 

In this non-perturbative regime, the description of hadrons is quite far 
removed from the description of the gluonic and quark fields in the QCD 
Lagrangian. Because only colour-singlet states survive, the hadrons are all 
composites of quarks and gluons. One example emphasises this: the nucleon 
has a mass which is very much greater than the sum of the quark masses 
of the three valence quarks comprising it. This extra mass comes from the 
gluonic interactions of QCD. Another way to view this is that the naive 
quark model is a useful phenomenological tool but has constituent quarks 
with masses much greater than the QCD masses (ie masses as defined in 
the Lagrangian). It is important to understand why this is approximately 
what QCD requires and to find where QCD departs from the naive quark 
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model. 

One way to characterise the manner in which QCD goes beyond the 
naive quark model is through the concept of exotic states. Here exotic is 
taken to mean 'not included in the naive quark model'. In order to discuss 
exotic states, we need to summarise what the naive quark model contains. 
Basically the degrees of freedom are the valence quarks (ie quark-antiquark 
for a meson and 3 quarks for a baryon) with masses and interactions given 
by some effective interaction. The consequences of this are that only certain 
J PC values will exist and that the number of states with different quark 
flavours is specified. So, concentrating on mesons made of the three flavours 
of light quarks (u, d, s), one expects a nonet of mesons with the flavours 
(ud, du, uu ± dd, ss, us, ds, su, sd). This is indeed what is found for 
vector mesons (p, u, <j), K*). It is also possible within the quark model for 
the flavour-singlet states (uu + dd, ss) to mix, as found for the pseudoscalar 
mesons. What would be exotic is for a tenth state to exist. For mesons 
with orbital angular momentum L between the quark and antiquark the 
allowed J PC values are shown below. Thus spin-parity combinations such 
as CP - ,!] -1 ,1 h , 2 H are termed spin-exotic since they cannot be made 
from a quark plus antiquark alone. 
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It has been a considerable challenge to build a machinery that allows 
non-perturbative calculations in QCD with all systematic errors deter- 
mined. The most controlled approach to non-perturbative QCD is via 
lattice techniques in which space-time is discretized and time is taken as Eu- 
clidean. The functional integral is then evaluated numerically using Monte 
Carlo techniques. 

Lattice QCD needs as input the quark masses and an overall scale 
(conventionally given by Aqcd)- Then any Green function can be eval- 
uated by taking an average of suitable combinations of the lattice fields 
in the vacuum samples. This allows masses to be studied easily and ma- 
trix elements (particularly those of weak or electromagnetic currents) can 
be extracted straightforwardly. Unlike experiment, lattice QCD can vary 
the quark masses and can also explore different boundary conditions and 
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sources. This allows a wide range of studies which can be used to diagnose 
the health of phenomenological models as well as casting light on experi- 
mental data. 

Here we will concentrate on lattice results for the spectrum of spin-exotic 
mesons: which are commonly called "hybrid mesons" since they must have 
additional degrees of freedom (such as gluonic) to achieve those spin-parity 
values. 

One limitation of the lattice approach to QCD is in exploring hadronic 
decays because the lattice, using Euclidean time, has important contribu- 
tions from low lying thresholds 1 which can obstruct the study of decay 
widths. The finite spatial size of the lattice implies that two-body states 
are actually discrete. By measuring their energy very precisely as the spa- 
tial volume is varied, it is possible 2 to extract the scattering phase shifts 
and hence decay properties. For on-shcll transitions, it is possible to es- 
timate hadronic transition strengths more directly 3 and this has recently 
been checked for the case of p meson decay to nn. This approach has been 
used to explore hybrid meson decay rates. 

2. Hybrid Mesons 

A hybrid meson is a meson in which the gluonic degrees of freedom are 
excited non-trivially. The most direct sign of this would be a spin-exotic 
meson, since that could not be created from a qq state with unexcited 
glue. A spin-exotic meson could, however, be a qqqq or meson-meson state 
and that possibility will be discussed. I first discuss hybrid mesons with 
static heavy quarks where the description can be thought of as an excited 
colour string. The situation concerning light quark hybrid mesons is then 
summarised 

2.1. Heavy quark hybrid mesons 

This topic has a long history: the first paper 4 was published 20 years ago. 
The basic ideas have not changed, and I summarise them here. 

Consider QQ states with static quarks in which the gluonic contribution 
may be excited. We classify the gluonic fields according to the symmetries 
of the system. This discussion is very similar to the description of elec- 
tron wave functions in diatomic molecules. The symmetries are (i) rotation 
around the separation axis z with representations labelled by J z (ii) CP 
with representations labelled by g(+) and u(— ) and (iii) C3?. Here C in- 
terchanges Q and Q, P is parity and 3? is a rotation of 180° about the 
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Figure 1. The potential energy between static quarks at separation R (in units of ro ~ 
0.5 fm) 6 . The symmetric gluonic field configuration is shown by the lower points while 
the IT U excited gluonic configuration is shown above. The energy levels in these potentials 
for b quarks are shown using the adiabatic approximation. 



mid-point around the y axis. The CSft operation is only relevant to classify 
states with J z = 0. The convention is to label states of J z = 0,1,2 by 
X, n, A respectively. The ground state will have J z — and CP = +. 

The exploration of the energy levels of other representations has a long 
history in lattice studies 4 ' 6 . The first excited state is found to be the 
This can be visualised as the symmetry of a string bowed out in the x 
direction minus the same deflection in the — x direction (plus another com- 
ponent of the two-dimensional representation with the transverse direction 
x replaced by y), corresponding to flux states from a lattice operator which 
is the difference of U-shaped paths from quark to antiquark of the form 

n - u. 

The picture of the gluon flux between the static quarks suggests that the 
excited states of this string may approximate the excited potentials found 
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from the lattice. In the simplest string theory, the first excited level has 
H u symmetry and is at energy ir/R above the ground state. This is indeed 
approximately valid and a closer approximation is to use a relativistic ver- 
sion 5 (namely E m (R) = {a 2 R 2 + 2ir<r(m - 1/12)) 1 / 2 for the m-th level), 
see also ref. 7 for a recent comparison of this expression. 

Recent lattice studies 7 have used an asymmetric space/time spacing 
which enables excited states to be determined comprehensively. These re- 
sults confirm the finding that the Ii u excitation is the lowest lying and 
hence of most relevance to spectroscopy. 

From the potential corresponding to these excited gluonic states, one can 
determine the spectrum of hybrid quarkonia using the Schrodingcr equation 
in the Born-Oppenhcimer approximation. This approximation will be good 
if the heavy quarks move very little in the time it takes for the potential 
between them to become established. More quantitatively, we require that 
the potential energy of gluonic excitation is much larger than the typical 
energy of orbital or radial excitation. This is indeed the case , especially 
for b quarks. Another nice feature of this approach is that the self energy of 
the static sources cancels in the energy difference between this hybrid state 
and the QQ states. Thus the lattice approach gives directly the excitation 
energy of each gluonic excitation. 

The n M symmetry state corresponds to excitations of the gluonic field in 
quarkonium called magnetic (with L PC = 1^ ) and pseudo-electric (with 
1 '") in contrast to the usual P-wave orbital excitation which has L PC = 
Thus we expect different quantum number assignments from those 
of the gluonic ground state. Indeed combining with the heavy quark spins, 
we get a degenerate set of 8 states: 
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Note that of these, J PC = 1 '", H and 2 H are spin-exotic and hence 
will not mix with QQ states. They thus form a very attractive goal for 
experimental searches for hybrid mesons. 

The eightfold degeneracy of the static approach will be broken by various 
corrections. As an example, one of the eight degenerate hybrid states is a 
pseudoscalar with the heavy quarks in a spin triplet. This has the same 
overall quantum numbers as the S-wave QQ state (t]b) which, however, has 
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the heavy quarks in a spin singlet. So any mixing between these states 
must be mediated by spin dependent interactions. These spin dependent 
interactions will be smaller for heavier quarks. It is of interest to establish 
the strength of these effects for b and c quarks. Another topic of interest 
is the splitting between the spin exotic hybrids which will come from the 
different energies of the magnetic and pseudo-electric gluonic excitations. 

One way to go beyond the static approach is to use the NRQCD approx- 
imation which then enables the spin dependent effects to be explored. One 
study 7 finds that the L PC = 1 H and 1 h excitations have no statistically 
significant splitting although the 1 H excitation does lie a little lighter. 
This would imply, after adding in heavy quark spin, that the J PC = 1 
hybrid was the lightest spin exotic. Also a relatively large spin splitting was 
found 8 among the triplet states considering, however, only magnetic glu- 
onic excitations. Another study 16 explores the mixing of non spin-exotic 
hybrids with regular quarkonium states via a spin-flip interaction using 
lattice NRQCD. 

Confirmation of the ordering of the spin exotic states also comes from 
lattice studies with propagating quarks 9 < 10 * 14 which are able to measure 
masses for all 8 states. I discuss that evidence in more detail below. 

Because of the similarity of the lightest hybrid wavefunction with that 
of the 2S state (which has L = 1), it is convenient to quote mass differences 
between these states. Within the quenched approximation, the lattice evi- 
dence for bb quarks points to a lightest hybrid spin exotic with J PC = 1 h 
at an energy given by (mjj — m 2 s)r =1.8 (static potential 6 ); 1.9 (static 
potential 7 , NRQCD 8 ); 2.0 (NRQCD 7 ). These results can be summarised 
as (ma — m 2 s)r = 1.9 ± 0.1. Using the experimental mass of the T(25), 
this implies that the lightest spin exotic hybrid is at mn = 10.73(7) GeV 
including a 10% scale error. Above this energy there will be many more 
hybrid states, many of which will be spin exotic. 

The results from a study with Nj = 2 flavours of sea-quarks show very 
little change in the static potential as shown by SESAM 11 and as illustrated 
in fig. 2 and also relatively little change in NRQCD determinations 8 of 
mass ratios such as (mjj — TO2s)/(toip — mis). Expressed in terms of tq 
(using ro = 1.18 /y/a) this gives (m# — m2s) r a = 2.4(2), however. This is 
significantly larger than the quenched result and, using the IP — IS mass 
difference to set the scale, yields a prediction 8 for the lightest hybrid mass 
of 11.02(18) GeV. 
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2.2. Hybrid meson decays 

Within this static quark framework, one can explore the decay mechanisms. 
One special feature is that the symmetries of the quark and colour fields 
about the static quarks must be preserved exactly in decay, hence the light 
quark-antiquark pair produced must respect these symmetries. This has the 
consequence that the decay from a IT M hybrid state to the open-6 mesons 
(BB, B*B, BB*, B*B*) will be forbidden 19 ' 21 if the light quarks in the 
B and B* mesons are in an S-wave relative to the heavy quark (since the 
final state will have the light quarks in either a triplet with the wrong CP 
or a singlet with the wrong J z where z is the interquark axis for the heavy 
quarks). The decay to i?**-mesons with light quarks in a P-wave is allowed 
by symmetry but not energetically. 
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Figure 2. The potential energy V(R) (in lattice units with a=0.0972 fm) versus quark 
separation R in fm for 2 flavours of sea quark. The energies are given for the ground 
state and first excited gluonic state and for the two body state of ground state potential 
plus scalar meson (/o) in a P-wave with the minimum non-zero momentum. The on-shcll 
transition can be evaluated when R as 0.2 fm. 
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In the heavy quark limit, the only allowed decays are when the hybrid 
state de-excites to a non-hybrid state with the emission of a light quark- 
antiquark pair. Since the U u hybrid state has the heavy quark-antiquark 
in a triplet P-wave state, the resulting non-hybrid state must also be in a 
triplet P-wave since the heavy quarks do not change their state in the limit 
of very heavy quarks. Thus the decay for b quarks will be to Xb + M where 
M is a light quark-antiquark meson in a flavour singlet. This proceeds 
by a disconnected light quark diagram and it would be expected 20 that 
the scalar or pseudoscalar meson channels are the most important (ic they 
have the largest relative OZI-rule violating contributions). This transition 
can be estimated on a lattice when the intial and final energies are similar. 
This is the case 21 for the II^ de-excitation to ground state gluonic field plus 
/o meson when the inter quark separation is around 0.2 fm which allows a 
lattice evaluation of the hadronic transition strength - see fig. 2. Indeed the 
dominant mode (with a width of around 100 MeV) is found to be with M as 
a scalar meson, namely H — > \b + fo, whereas when M is an r\ or rf meson 
the transition strength is less that a few MeV. There will be modifications 
to this analysis coming from corrections to the heavy quark limit (of order 
1 / niQ where tuq is the heavy quark mass) which might allow hybrid meson 
transitions to BB, etc, but these have not been evaluated yet. 

In this heavy quark (or static) limit, the spin-exotic and non spin-exotic 
hybrid mesons are degenerate. For the latter, however, the interpretation 
of any observed states is less clear cut, since they could be conventional 
quark antiquark states. Moreover, the non spin-exotic hybrid mesons can 
mix directly (ic without emission of any meson M) with conventional quark 
antiquark states once one takes into account corrections (of order 1/Mq) to 
the static approximation applicable for heavy quarks with physical masses. 

It is encouraging that the decay width comes out as relatively small, 
so that the spin-exotic hybrid states should show up experimentally as 
sufficiently narrow resonances to be detectable. This decay analysis does 
not take into account heavy quark motion or spin-flip and these effects will 
be significantly more important for charm quarks than for 6-quarks. 

2.3. Light quark hybrid mesons 

I now focus on lattice results for hybrid mesons made from light quarks 
using fully relativistic propagating quarks. There will be no mixing with 
qq mesons for spin-exotic hybrid mesons and these are of special interest. 
The first study of this area was by the UKQCD Collaboration 9 who used 
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operators motivated by the heavy quark studies referred to above to study 
all 8 J values coming from L PC = l" 1 and 1 h excitations. The result- 
ing mass spectrum gives the J PC = 1 h state as the lightest spin-exotic 
state. Taking account of the systematic scale errors in the lattice deter- 
mination, a mass of 2000(200) MeV is quoted for this hybrid meson with 
ss light quarks. Although not directly measured, the corresponding light 
quark hybrid meson would be expected to be around 120 MeV lighter. 

A second lattice group has also evaluated hybrid meson spectra with 
propagating quarks from quenched lattices. They obtain 10 masses of the 
1 h state with statistical and various systematic errors of 1970(90)(300) 
MeV, 2170(80) (100) (100) MeV and 4390(80) (200) MeV for nn, ss and cc 
quarks respectively. For the H spin-exotic state they have a noisier signal 
but evidence that it is heavier. They also explore mixing matrix elements 
between spin-exotic hybrid states and 4 quark operators. 

The first analysis 14 to determine the hybrid meson spectrum using full 
QCD used Wilson quarks. The sea quarks used had several different masses 
and an extrapolation was made to the limit of physical sea quark masses, 
yielding a mass of 1.9(2) GcV for the lightest spin-exotic hybrid meson, 
which again was found to be the 1 h . In principle this calculation should 
take account of sea quark effects such as the mixing between such a hybrid 
meson and qqqq states such as rjir, although it is possible that the sea quark 
masses used are not light enough to explore these features. 

A recent dynamical quark study from 2+1 flavours of improved stag- 
gered quarks has also produced results 12 . They also compare their results 
with quenched calculations and find no significant difference, except that 
the ambiguity in fixing the lattice energy scale is better controlled in the 
dynamical simulation since different reference observables are closer to ex- 
periment. Their summary result for the 1 h hybrid with strange quarks 
is 2100 ± 120 MeV, in agreement with earlier results. They note that the 
energies of two-meson states (such as n + bi or K + K(l + ) ) with the hybrid 
meson quantum numbers are close to the energies they obtain. This sug- 
gests that these two-particle states, which are allowed to mix in a dynamical 
quark treatment, may be influencing the masses determined. A study of 
hybrid meson transitions to two particle states is needed to illuminate this 
area, using techniques such as those used for heavy quark hybrid decay 21 
and decays of light quark vector mesons 3 . 

The lattice calculations 940,144245 Q f ^ e light hybrid spectrum arc 
in good agreement with each other. They imply that the natural energy 
range for spin-exotic hybrid mesons is around 1.9 GeV. The J PC = 1~ + 
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state is found to be lightest. It is not easy to reconcile these lattice results 
with experimental indications 17 for resonances at 1.4 GeV and 1.6 GeV, 
especially the lower mass value. Mixing with qqqq states such as rin is 
not included for realistic quark masses in the lattice calculations. Such 
effects of pion loops (both real and virtual) have been estimated in chiral 
perturbation theory based models 18 and they could potentially reconcile 
some of the discrepancy between lattice mass estimates (with light quarks 
which are too heavy) and those from experiment. This can be interpreted, 
dependent on one's viewpoint, as either that the lattice calculations are 
incomplete or as an indication that the experimental states may have an 
important meson-meson component in them. 

The light quark technique of using relativistic propagating quarks can 
also be extended to charm quarks, as was note above 10 . Another group has 
explored the charm quark hybrid states also using a fully relativistic action, 
albeit with an anisotropic lattice formulation 13 . Their quenched study is 
in agreement with the isotropic lattice result quoted above, finding a mass 
value of 4.428(41) GeV in the continuum limit for the 1 h hybrid where 
the scale is set by the 1 P\ — IS mass splitting (458.2 MeV experimentally) 
in charmonium. Their result is also consistent with that from NRQCD 
methods 8 applied to this case. These results all have the usual caveat 
that in quenched evaluations the overall mass scale of the energy difference 
from the IS* state at 3.067 GeV is uncertain to 10% or so (for example the 
(2S — 1S)/( 1 Pi — IS) is found to be 15% higher than experiment) which is 
a major source of systematic error (approximately ±140 MeV). They also 
produce estimates for other charmonium spin-exotic states: 0^ at 4.70(17) 
GeV and 2+" at 4.89(9) GeV. The 0"" state is not resolved. 

Thus masses near 4.4GeV are found for the charmonium 1 h state using 
relativistic quarks. The non-relativistic approach using NRQCD is expected 
to have big systematic errors for quarks as light as charm, but results 8 do 
agree with this value. The heavy quark effective theory approach has a 
leading term which corresponds to a static heavy quark, resulting in an 
estimate 6 of the spin-exotic charm state mass of 4.0 GeV. Here again the 
systematic error is potentially large for charm quarks. 

3. Conclusions and Outlook 

For hybrid mesons, there will be no mixing with qq for spin-exotic states 
and these are the most useful predictions. The J PC = 1~ + state is expected 
in the range 10.7 to 11.0 GeV for b quarks, 2.0(2) GeV for s quarks and 
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1.9(2) GeV for u, d quarks. Mixing of spin-exotic hybrids with qqqq or 
equivalently with meson-meson is allowed and will modify the predictions 
from the quenched approximation. A first lattice study has been made of 
hybrid meson decays. For heavy quarks, the dominant mode is string de- 
excitation to x + fo where f is a flavour singlet scalar meson (or possibly 
two pions in this state) . The magnitude of the decay rate is found to be of 
order 100 MeV, so this decay mode should still leave a detectably narrow 
resonance to be observed. It will be possible to explore hadronic decays for 
light-quark hybrid mesons in the near future, and this may help to elucidate 
the nature of the experimental spin-exotic signals found at 1.4 and 1.6 GeV. 
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